1 



o 
o 

(N 
o 

Q 



Progress on R Measurement through ISR with BABAR 



Michel Davier 

Laboratoire de l'Accelerateur Lincaire 
IN2P3/CNRS ct Universite de Paris-Sud 
91898 Orsay, France 
E-mail: davier@lal.in2p3.fr 



> 

o 

m 
O 

">< 

<D 
Oh! 



X 



The status of the ISR analysis with BABAR is presented. Preliminary results are presented on the process 
e + e _ — * 27r + 27r~ and J/V> decays including a determination of its total width and the branching ratios into 
27T+27T", K + K~TY + n~ and 2K + 2K~ modes. 



1. Introduction 

High- luminosity e + e~ storage rings designed 
for the study of CP violation, such as KEK-B, 
PEP-B, and DA$NE, have opened a new possi- 
bility to measure e + e~ annihilation using initial- 
state radiation (ISR) flTM\. A large mass range 
becomes accessible in a single experiment, espe- 
cially at B factories, contrary to the case with 
fixed energy machines which are optimized only 
in a limited region. The broad-band coverage may 
result into a better understanding of systematic 
effects since only one experimental set-up is in- 
volved. 

Apart from their own physics interest, rang- 
ing from hadronic spectroscopy to tests of QCD, 
measurements of the ratio R = a(e + e~ — > 
hadrons)/<7 p t provide the necessary input to dis- 
persion integrals for computing hadronic vacuum 
polarization contributions. This is indeed the 
case when calculating the theoretical prediction 
for the muon magnetic anomaly and the run- 
ning of the electromagnetic coupling to the Z 
mass, a(Af|), one of several ingredients in pre- 
cision tests of the electroweak theory bearing on 
such fundamental issues as the structure of the 
Higgs sector. 

The cross section for the ISR process e + e~ — > 
/7 with a particular final state / is related to 
the cross section <r/(s) for the direct annihilation 
e + e~ — > / through 



da(s, x) 
dx 



where x — 2E* / yfs, E* being the radiated pho- 
ton energy in the e + e _ frame and yfs the nominal 
centre-of-mass energy of the collider. The quan- 
tity s' = s(l — x) represents the mass squared of 
the final-state / system. The radiator function 
W(s, x) providing the virtual photon spectrum 
can be computed including radiative corrections 
to an accuracy better than 1% 0|. The produc- 
tion of ISR photons is strongly peaked along the 
initial beams, but for yfs ~ 10 GeV the fraction 
at large angle is relatively large, i.e. 10 — 15% in 
the BABAR central detector. 

The measurement of the corresponding leptonic 
process e + e~ — * /i + /i~7 provides the ISR lumi- 
nosity. Thus the Born cross section crj(s') is ob- 
tained through 
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where ANf^ is the number of detected fj events 
in the bin of width As' centred at s', e/ is the 
detection efficiency for the final state / and Sp SR 
is the fraction of the cross section for e + e~ — > / 7 
where the photon is emitted by the final-state 
particles. The latter quantity can be sizeable 
for the mi channel, but negligible for most of 
the low-energy hadronic states which have van- 
ishingly small cross sections at the nominal en- 
ergy yfs. The radiative corrections for the initial 
state and vacuum polarization effects cancel in 
the ratio J2J, as does the photon efficiency. 
The BABAR detector at the PEP-II asymmet- 



2 



ric e + e~ storage ring is well suited for this study. 
The different triggers routinely used allow one to 
select ISR processes with a 95 — 98% efficiency. A 
complete description of BABAR can be found else- 
where [5] . All main detector components are used 
in this analysis as particle identification plays a 
crucial role in separating the different channels, 
ft is also important in correctly computing the 
s' value for each event as it is obtained from the 
produced particles momenta since the photon en- 
ergy is rather insensitive to s' in the low-energy 
regime of interest. In particular, K-7T separation 
in quartz radiators (DIRC) is an essential ingre- 
dient. The information from the tracking sys- 
tem (silicon vertex detector and drift chamber) is 
used to measure the angles and the momentum 
of the charged particles. Muons are identified in 
the instrumented flux return (IFR) of the super- 
conducting solenoidal magnet, while photons are 
detected in the fine-grain Csl calorimeter. 

The data used in this analysis correspond to an 
integrated luminosity of 89.3 fb _1 , collected both 
at the T(4S') and in the nearby continuum. 



2. Event preselection 

ISR events are preselected requiring a large- 
angle photon in the central part of the detector 
with an energy in the centre-of-mass larger than 
3 GeV. This requirement is loosened down to 0.5 
GeV for 2-track events which have at least one 
track identified as a muon, in order to collect a 
large statistics of muon tracks for calibrating the 
muon-idcntification performance. The number of 
tracks pointing well to the interaction region is re- 
quired to be even and larger than two, and their 
total charge should be equal to zero. A loose kine- 
matic constraint is applied by asking the ISR pho- 
ton candidate to be inside a 0.3 radian cone whose 
axis is along the missing momentum vector con- 
structed from all the recorded recoiling particles. 

Essentially all low-energy exclusive ISR final 
states are retained by these general selection cri- 
teria, while about 2% of all BABAR events are 
kept. 



3. The di-muon final state 

The process e + e~ — > /i + /i~7 plays an impor- 
tant role in this analysis. For /i/i masses less 
than a few GeV, it provides the measurement of 
the ISR luminosity, while for masses larger than 
2 GeV, where hadron pairs are negligible, it en- 
ables one to measure the muon identification effi- 
ciency. For this latter purpose, a muon candidate 
is tagged by the opposite identified muon. The 
large statistics of muons collected in this way ef- 
ficiently covers the acceptance of the detector so 
that fine-grained efficiency maps in momentum, 
polar and azimutal angles can be constructed. 

A 1C kinematic fit using the two muon 
tracks effectively removes background from t + t~ 
events. After the mass and x 2 cuts, the pu- 
rity of the muon sample is found to be better than 
0.999. In addition, the fit improves the fifi mass 
resolution which decreases from 16 MeV down to 
8 MeV at the mass. 

Fig. n shows the derived spectrum of ISR lu- 
minosity per 100 MeV energy bins. Thus, the 
BABAR luminosity of 89.3 fb^ 1 corresponds to 
an e + e~ energy scan taken all together with 0.7 
pb-yi00 MeV at 1 GeV and 3.6 pb-yi00 MeV 
at 4 GeV. It provides a statistically competitive 
hadron sample, especially in the range between 
1.4 and 3 GeV. 

4. Hadronic channels 

The main approach to the R measurement is 
the study of exclusive hadronic channels, con- 
strained by a kinematic fit and with identified 
charged particles. Work is in progress along these 
lines, including the measurement of the follow- 
ing final states: 7r + 7r _ , K + K~ , pp, ir + ir~ir°, Ait, 
57T, 6-7T, 7T7r?7, KK?t, KKnir, 2K2K, KKr), with 
charged and neutral (K®) kaons. Another ap- 
proach, based on an inclusive method where the 
hadronic mass is deduced from the ISR photon 
energy, is also in progress, with the difficulty that 
the resolution deteriorates rapidly for low recoil 
masses. 

The study of the ir + ir~ final state is an impor- 
tant part of the programme, especially in view 
of persisting discrepancies between existing e + e~ 
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Figure 1. The differential ISR luminosity in 
100-MeV bins obtained from the measurement of 
e + e~ — > (i + /i~"f. The points near the p(770) 
and J/0 resonances are excluded, because of pion 
misidentiftcation as muon in the first case and the 
non-QED nature of the second process. 



results and isospin-corrected r data [61 7j . Such a 
study requires a good control of systematic un- 
certainties at the 1% level or better. The ISR 
method has several advantages compared to the 
standard fixed-energy annihilation: (i) a large 
mass range is simultaneously covered in a single 
experiment, (ii) the cross section can be measured 
down to threshold with excellent acceptance, and 
(iii) easier particle identification can be achieved, 
thanks to the Lorentz boost given to the hadronic 
system. 

Since the 7r + 7r~ and final states have the 

same topology, it is advantageous to directly mea- 
sure the ratio N^/N^ since many effects cancel: 
the e + e~ luminosity, ISR and vacuum polariza- 
tion radiative corrections, the photon detection 
efficiency. Small corrections still need to be ap- 
plied concerning trigger efficiencies (calorimeter 
triggers are more efficient for pions), tracking ef- 
ficiency (because pions can produce secondary in- 



teractions in the tracking detector), FSR correc- 
tions. All these effects can in fact be studied di- 
rectly on the data. The major experimental task 
remains the determination of the particle identi- 
fication efficiency matrix, using pure samples of 
known particles, as outlined above for the muons, 
as a function of momentum and the 2-dimensional 
location of the particle in the relevant detectors. 
Such a detailed study is still in progress. 

Preliminary results are available for the fi- 
nal state with four charged particles, namely 
27T+27T-, K+K-tt+it-, and 2K+2K-. Fig. 
shows the mass distribution of the events satis- 
fying the 27r + 27r - kinematics and K/7r identifica- 
tion constraints. No muon identification was ap- 
plied as no corresponding background is expected 
in this four-track sample. The data are seen to be 
dominated by a broad peak at 1.5 GeV originat- 
ing from the p(1450) and p(I700) resonances, a 
shoulder near 2 GeV, and sharp peaks at the J/ip 
and ip' masses. The latter signal is the result of 
the decay chain 0' — * 7r + 7r~ J/V>, 3/ip — > 
The background from other hadronic channels is 
small and subtracted in each mass bin. Overall, 
about 70k events are obtained leading to small 
statistical uncertainties. The normalization un- 
certainty is estimated to be 5%, dominated by 
the luminosity determination and the dependence 
of the detection efficiency on the dynamics (be- 
low 2 GeV the intermediate state it ai is found 
to be dominant). The derived cross section is 
given in Fig. [3] The accuracy of the BABAR data 
is comparable to the latest most precise results 
from CMD-2 8 and SND JH] below 1.4 GeV and 
the agreement is good. Between 1.4 and 2 GeV 
the quality of the BABAR data exceeds that from 
DCI and Adone. The range above 2 GeV is only 
covered by BABAR. 

One can illustrate the impact of the BABAR re- 
sults on the 2tt + 2it~ annihilation cross section 
obtained through ISR by computing its contribu- 
tions to the lowest-order hadronic component of 
the muon magnetic anomaly. Using the same en- 
ergy range from threshold up to 1.8 GeV as in 
Ref. |l()j for a convenient comparison, one gets 

a had,LO [47r] = (i2.95±0.64 eiEp ±0.13 ra( i) 1(T 10 (3) 
to be compared to the values JU| obtained us- 
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Figure 2. The J^-pion invariant mass distribu- 
tion for kinematically constrained events of the 
type e + e~ — > 2it + 2'k~j. No muon/pion separa- 
tion is applied. The points with error bars indi- 
cate the level of background obtained from a study 
of the x 2 distribution of the kinematical fit. The 
cross-hatched histogram corresponds to the non- 
ISR background as given by the JETSET simu- 
lation of annihilation events. Sharp signals are 
seen at the J/ip and ip' masses. 



ing all previously available e + e data, (14.21 ± 



0.87 exp ± 0.23, 
difficult t~ — 



2 d) 10 , and from r data on the 
7r _ 37r°^ r decay mode, (12.35 ± 



0.96 exp ± 0.40 sc/(2) ) 10~ 10 . In the above results, 
the uncertainties have been separated into exper- 
imental sources 'exp' (statistics and systematics), 
radiative corrections 'rod' and isospin-breaking 
corrections 'SU{2)\ The BABAR results lead to 
a significant improvement in this channel. 

5. Radiative return to the J/tp 

For a narrow state of mass M such as J/ip, 
decaying into a final state /, one has 



12n 2 T ee B f 
°"J/V = ^f s W{s,xo) 



(4) 



Figure 3. The e + e — > 2tt + 2tt cross section 
obtained from ISR at BABAR in comparison with 



all existing 



data. 



with x — l — M 2 /s. The total J/ip ISR cross sec- 
tion amounts to 0.036 nb for s = (10. 58) 2 GeV 2 , 
leading to 4 10 6 events for the present 130 fb _1 
of BABAR luminosity (the angular acceptance of 
the detected ISR photon is about 0.1). 

For the final state, one can directly com- 

pare the rate in the 3/tjj peak (about 8000 events) 
to the continuum (see Fig.0J) and derive the value 
of the product 



(0.330 ± 0.008 stat ± 0.007 svst ) keV 



(5) 



The systematic error receives contributions from 
the background estimation, the J/ip line shape, 
the radiative corrections and Monte Carlo statis- 
tics. Using the world averages for Bj^^^+ |U - 
and Bj /^_, e + e - , one can derive the J/i/> electronic 
and total widths, T 3/ ^ e+e - = (5.61 ± 0.20) keV 
and Tj/^ = (94.7 ± 4.4) keV, to be compared to 
the world average values of (5.26 ± 0.37) keV 
and (87 ± 5) keV, respectively. 
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Figure 4. The (J, + (i invariant mass distribution 
in the vicinity of the J/ip resonance. 
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From the mass spectrum given in Fig. [21 the 
following quantities can be derived: 



(1.95±0.14 stat ±0.13 a „ st 

--Bj/i/'->M + M~-^V''- > J/' ! /"r+ir 



■47T 



10" 2 keV 



L ip' — >-e+ e 



(6) 
(7) 



(4.50 ± 0.18 stat ± 0.22 syst ) 10- 2 keV 

Using the world averages JT] for the 3/tp and ?// 
electronic widths, the results Bj/^^ = (3.70 ± 
0.27 ± 0.36) 10~ 3 and B^,^. 3 /^ +ir - = 0.361 ± 
0.015 ±0.037 are obtained, in agreement with the 
world average values JT] of (4.0 ± 1.0) 10 -3 and 
0.305 ± 0.016, respectively. 

Taking advantage of the K/ir Cherenkov iden- 
tification, one can measure the J/ip decay rate 
into final states involving K + K~ pairs. The 
corresponding mass spectra are given in Figs. 
and El Using again the world average value 
for the electronic width, the following branch- 
ing ratios are measured: Bj/^/f+jc-i+i- = 
(6.25±0.50±0.62) 10~ 3 and B 3/4 ^ K + K - K+K - = 
(6.9 ± 1.2 ± 1.1) 10~ 4 . Both determinations are 
more precise than the current world averages |1 lj 
of (7.2±2.3) 10~ 3 and (7.±3.) 10~ 4 , respectively. 



Figure 5. The K + K~tt + it~ mass distribution for 
ISR-produced e + e~ — > K + K~tt + tt~ in the J/ij) 
region. 
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Figure 6. The K + K~ K + K~ mass distribution 
for ISR-produced e + e~ — > K + K~ K + K~ in the 
J/i/j region. 
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6. Conclusions 

Preliminary results from BABAR on final states 
produced through ISR show the high physics po- 
tential of this sample which should yield pre- 
cise measurements of e + e~ annihilation cross sec- 
tions. The ISR method will allow one to cross 
check the results of the 'direct' experiments, with 
the significant advantage that the important en- 
ergy range from threshold to 3-4 GeV can be cov- 
ered in a single experiment, thereby reducing sys- 
tematic uncertainties. 

The ratio R will be measured from the sum 
of exclusive channels, providing input for vacuum 
polarization calculations, QCD studies and inves- 
tigations in hadronic spectroscopy. Amazingly, 
the current data is already providing very com- 
petitive results for i/ip and ip' decays. 
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